(N 



% 



Nonparametric estimation of multivariate extreme-value 

copulas 



a,l,* 



Gordon Gudendorf '\ Johan Segers 

f^*) ' "Institut de statistique, biostatistique et sciences actuarielles, Voie du Roman Pays 20, B-1348 

O^l . Louvain-la-Neuve, Belgium. 

> 

o 

:^ 

0> ; Abstract 

Extreme-value copulas arise in the asymptotic theory for componentwise maxima 
\J^ • of independent random samples. An extreme- value copula is determined by its 
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certain shape constraints that arise from an integral transform of an underlying 
measure called spectral measure. Multivariate extensions are provided of certain 
rank-based nonparametric estimators of the Pickands dependence function. The 
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^ ■ is enforced by replacing an initial estimator by its best least-squares approxima- 

"^ . tion in the set of Pickands dependence functions having a discrete spectral mea- 
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CN \ estimators is demonstrated and the finite-sample performance of the estimators is 

r^ I investigated by means of a simulation experiment. 
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1. Introduction 

Extreme-value copulas arise in the asymptotic theory for componentwise max- 
ima of independent random samples. They provide the dependence structures for 
the class of multivariate extreme-value or max-stable distributions. More gener- 
all y, they constitute a flexible cla ss of models for describing positive association; 
see lGudendorf and SegersI (120101) for a survey. 

In this paper we will focus on the nonparametric estimation of extreme- value 
copulas in general dimensio ns. In particular, we aim a t multivariate extensions of 



the rank -based estimators inlGenest and SegersI (120091) and the projection method- 



ology in lFils-Villetardetal.1 (120081) . 

Let Xj = {Xji, . . . ,Xj p), i e {1, ... ,n}, be an independent random sample 
from a 7?-variate, continuous distribution function F with margins Fi, . . . ,Fp and 
copula C, that is, 

F{x) = CiF,(xi),...,Fp(Xp)l xeW, 

where F(x) = P(X ^ x) (componentwise inequalities), Fj(xj) = F(Xj ^ Xj), 
and C is the joint distribution function of (Fi(Xi), . . . , Fp(Xp)). We are interested 
in nonparametric estimation of C in the model where the margins Fi,. . . ,Fp are 
completely unknown (but continuous) and C is known to be an extreme-value 
copula. 

A p-variate copula C is an extreme-value copula if there exists a finite Borel 

measure H on the unit simplex A^^i = {(vvi, . . . , Wp) s [0, 1]'^ : wi H 1- Wp = 1}, 

called spectral measure, such that 

C(h) = exp(-^(- log Ml,..., -log Mp)), Me(0, 1]'', (1.1) 

the tail dependence function t : [0, ooy -^ [0, oo) being given by 

t{x)= \ va2i^{v^x\,...,VpXp]H{diV). (1.2) 

The spectral measure H is arbitrary except for the p moment constraints 

vy//(dv) = l, 7 6(1,...,;.}, (1.3) 



A„ 



which are equivalent to requiring that the margins of C be uniform on (0, 1). 

The tail dependence function t in (11.21 ) is convex, homogeneous of order one, 
that is t(cx) = c £(x) for c > 0, and satisfies max(xi ,.. .,Xp) ^ £(x) ^ Xi + 



+ Xp for all jc 



G [0, oo)P. By homogeneity, i is characterized by the Pickands 
dependence function A : Ap_i — > [1/p, 1], which is simply the restriction of £ to 
the unit simplex: for x e [0, oo)p \ {0}, 



i{x) = (jci + • • • + Xp)A{wi, . . . , Wp-i) 



where w; = 



Xi + • • • + Xp 

Here and further on, we frequently identify Ap_i with {(wi, . . . , Wp-i) e [0, 1]''"' : 

Wi H h Wp_i ^1}. The extreme- value copula C can be expressed in terms of A 

via 

log Ml logMp_i 



C(u) = exp|(2y=i log My) A 
forw e(0, 1F\{(1,...,1)}, with 



25'=il0gM/""Zy=il0gMy 



Aiw) 



-I 

J A,, -I 



maxjviwi, . . . , VpWp}H{dv) 



W G Ap_i, 



(1.4) 



(1.5) 



see 



Pickandsl (|1981|) and IZhang et alJ (|2008h . The function A is convex as well 



and satisfies max(wi, . . . , w,,) ^ A(w) ^ 1 for all w e A„_ 



p-i- 



Nonparametric estimator s for A have initi ally been developed i n iPickands 



in 



(1198 11). with modificatio ns in lDeheuvelsl(ll99lh and lHall and Tajvidil (l2000l) . and 



Caperaa et alJ (| 19971 ). These estimators will be referred to as the Pickands 
and CFG estimators, respectively; see Section [3] for definitions. In the previ- 
ously cited papers, the marginal distributions were assumed to be known. The 
m ore realistic case of unknown margins has been treated in the bivariate case 
in I Jimenez et al.l (|200lb for a submodel and in lGenest and SegersI (|2009[) for the 



gen eral model. Multivariate exten sions have been proposed in lzhang et al.l (|2008h 
and lOudendorf and SegersI (|201l|) for the case of known margins. In Section[3l we 
will provide a proof for the convergence of these estimators in case of unknown 
margins being estimat ed by the empirical distrib ution functions, thus generaliz- 
ing the main results in iGenest and SegersI (l2009h to arbitrary dimensions. As in 
Kojadinovic and YanI (|2010h and lOenest et al.l (1201 lb . the estimators could also be 
used as a starting point for goodness-of-fit tests, but for brevity we do not pursue 
this here. Finally, a new type of nonparametric estimator has been proposed in 
Bucher et al.l (|201l|) for the bivariate case. 

In the proofs of the asymptotic normality of the Pickands and CFG est ima- 
tors, a certain expansion of the empirical copula process due to lStutd (|l984r) and 



Tsukaharal (l2005h plays a crucial role. The second-order derivatives of extreme- 
value copulas typically exhibit explosive behaviour near the cor ners of the hyp er- 
cube, violating the assumptions in the two papers just cited. In lSegersl(|201lb . it 
was shown that the same expansion continues to hold under much weaker con- 
ditions on the partial derivatives. In Section |2l these issues are considered for 
multivariate extreme- value copulas. 

As the estimators for A considered here fail to be Pickands dependence func- 
tions themselves, it is natural to ask how to enforce the shape constraints on such 
functions in the estimation procedure. In dimension p = 2,it is sufficient to ensure 
that A is convex and takes values i n the range maxjw ) ^ A(w) ^ 1, for instance 
by truncation and convexification (IDeheuvelsl 1 1 99 ih . In dimensi on p ^ 3, how- 
ever, this procedure is no longer sufficient (IBeirlant et al.L l2004l page 257) and 
one needs to rely on the spectral representation in (11.51). I n Section |4] we will 
apply an projection methodology (IFils-Villetard et al.l . |2008|) to obtain valid esti- 
mates: an initial estimate is replaced by its best least-squares approximation in the 
set of Pickands dependence functions corresponding to discrete spectral measures 
supported on a fine grid. 

The results of a simulation experiment aimed at investigating the finite-sample 
performance of the original and projected Pickands and CFG estimators are re- 
ported in Section [51 All proofs are relegated to the Appendix. 

Throughout the article, we will apply the following notations. For a space 
W, let £°°(W) and 'tfiW) denote the spaces of real-valued bounded and real- 
valued continuous functions respectively, where we endow both spaces with the 
uniform norm || • lU : / •-> sup^^;^- \f{x)\- Furthermore, the indicator function of 
the event E is denoted by 1(E). The arrow '-w' will stand for weak convergence. 
For any /j-variate real-valued function / with values in M, the first and second- 
order partial derivatives will be denoted by fi(x) = -^/(xi, . . .,Xp) and ftjix) = 



dxjdxj 



f(Xi, . . . ,Xp). 



2. Empirical copula processes 



Let Xi , X2, ... be an iid sequence of random vectors from a /)-variate multivari- 
ate distribution F with continuous margins Fi, . . ., Fp. If the margins Fi,. . . ,Fp 
are known, we can define the empirical cumulative distribution function C„ of the 
(unobservable) random sample i/,- = (t/,,i, . . . , t/,-p) = (Fi(X,j), . . . , Fp{Xi^p)) for 



z e {l,...,n}by 

1 " 
Cn(u) = -J]l(Ui,i^uu...,Ui,p^Up), He[0,lf, (2.1) 



with associated empirical process 

On = n^^^ (Cn - C) . 

For ease of notation, we will write 

anjiuj) = anil, ...,\,Uj,l,...,l) for J e {1, ...,p}. 



(2.2) 



(2.3) 



In practice, the marginal distributions will need to be estimated. If we are not 
willing to make any assumptions (except for continuity) we can estimate them by 
the empirical distribution functions 



1 " 
Fnj(x) = y liXij ^x), j £[!,..., p}, 

1=1 



(2.4) 



where we divided by n + 1 in order to avoid later problems at the borders. By so 
doing, we can construct n vectors f/,- = (t/;,i, . . . , t/,-,p) via 



U: 



I=\ 

for i e { 1 , . . . , n} and j e { 1 , . . . , />}. The empirical copula will be denoted by 

1 " 
Cn{u) = -Y^\[Ui,,i^uu...,Ui,p^Up), He[0,lF (2.6) 



(2.7) 



with associated empirical copula process 

C„ = n'l\Cn - C). 

In IStutd 119841) and iTsukaharal (|2005|) it was established that if all second- 
order derivatives of C exist and are continuous on [0, IY\ the processes a,, in (12. 2D 
and C„ in (|2.7I) are related via 

p 
C„(m) = «„(«) - X Cj{u) anjiuj) + Rniu), (2.8) 



the remainder term R„ satisfying 

m^ \Rn{u)\ = 0{n-^'\\ognyi^(\og\ognyi^) almost surely. (2.9) 

He[0,l]P 

Let i°°{\Q, l]P) denote the space of bounded real-valued functions on [0, 1]^, 
equipped with the topology of uniform convergence. Weak convergence of maps 



taking values in i°°{{0, 1]'') will be understood as in Ivan der Vaart and Wellnei 



(119961 Definition 1.3.3) and will be denoted by '-w'. By classical empirical pro- 
cess theory, we have a,, --* a as n ^ cx), the limiting process being a mean-zero 
tight Gaussian process with continuous trajectories and covariance function 

cov(Qr(H), a{v)) = C(u A v) - C{u) C(v), u,v e [0, 1]^, (2.10) 

where (u A v)j = mm(uj, Vj). In view of the expansion (|2.8I) . it then follows that in 
i°°([0, \]P), we have C„ -^ C as n ^ oo, where 

p 
C(m) = aiu) - ^ Cj(u) ajiuj) (2.1 1) 

i=i 

and cnjiuj) = a(l, . . . ,\,Uj,\, . . . , 1). 

Like many other copulas, extreme- value copulas do in general not have uni- 
formly bounded second-order partial derivatives. For instance, in the bivariate 
case, every copula having a positive coefficient of upper tail dependence will have 
first-order partial derivatives that fa il to have a continuous extension to the up- 



per corner (1, 1); see iSegersI (|201lL Example 1.1). As a consequence, the only 
bivariate extreme-value copula whose density is uniformly bounded is the inde- 
pendence copula. However, as shown in the same paper, for copulas satisfying 
Assumption [T] below, the expansion (I2.8I) -( |2.9I ) of the empirical copula process 
remains valid. 

Assumption 1. (CI) For every j e {\,...,p}, the first-order partial derivative 
Cj exists and is continuous on the set Vpj = {w 6 [0, 1]'' : < Uj < 1}. 

(C2) For every i,j£ {\, . . . ,p} (i and j not necessarily distinct), Cij exists and is 
continuous on Vpj n Vpj and 

sup max{M;(l - Ui),Uj(l - Uj)] \Cij{u)\ < oo. 

ueVpjriVpj 



In fact, for weak convergence C„ -w C in ^°°([0, 1]^) to hold, condition (CI) 
is already sufficient. In the context of multivariate extreme-value copulas, it will 
be of interest to have a readily verifiable condition on the stable tail dependence 
function € for Assumption[T]to hold. 

Assumption 2. (LI) For every j e {\, . . . ,p}, the first-order partial derivative -tj 
exists and is continuous on the set Wp^ = {x € [0, ooy : Xj > 0}. 

(L2) For every i, j £ {\, . . . ,p} (i and j not necessarily distinct), tij exists and is 
continuous on Wpj D Wpj and 

sup max(xi, xj) \'£ij(x)\ < oo. 

xeWpjnWpj 

Xl+---+Xp = l 

Proposition 1. For p-variate extreme-value copulas, (LI) implies (CI). If in ad- 
dition (L2) holds, then (C2) holds as well. 

In the bivariate case, sufficient conditions for (CI) and (C2) can be given in 
terms of the Pickands dependence function A(w) = £{\ - w, w), where w 6 [0, 1]: 
(CI) holds as soon as A is continuously differentiable on (0, 1), and (C1)-(C2) 
hold as soon as A is twice con tinuously differentiable on (0, 1) and supQ^,^,^j w(l - 



w)A"{w) < oo (lSegersl . l201lL Example 5.3). 

For completeness, we want to mention that in the above references, the empir- 
ical copula is not defined as in (12.61) but rather as 

Cf («) = F„(F^i(mi), . . . , F^p{up)), 

with F^Xuj) = mf{xj e M : F„j(xj) ^ uj}. Straightforward calculus shows that, 
in the absence of ties, 

sup |Cf (h) - Cn(u)\ ^ , 

He[0,l]'' n 



As a consequence, Stute's expansion (12.81) is valid for C„ if and only if it is valid 
for Cf . 



3. Nonparametric estimation of the dependence function 



Among th e most popular nonparametric estimators for A figur e the Pickands 



estima tor Al (|Pickandsl . ll98ll) and the estimator A^^'^ proposed by ICaperaa et al. 
(Il997h . referred to as the CFG estimator from now on. Writing 



^;(H') = A -. 



7=1 ' 

for w e Ap^i, with Uij as in (|2.5I) . these estimators are defined as 
1 



A^iw) 



- y Uw) and log Af ^Civ) = -- V logl.-dv) - y, 

1=1 1=1 



with 7 = 0.5772 . . . the Euler-Mascheroni constant. Explanations on the construc- 
tion of th ese estimators are provide d for instance in the or iginal references given 
before , in lGenest and SegersI (|2009r ) and in the survey paper lOudendorf and Segers 
( 20101) . The multivariate extension of the CFG estimator was presented in lzhang et al.l 
(|2008b . albeit under a different but equivalent form. 

In order to improve the small-sample properties of the above estimators, the 
endpoint constraints A(ej) = 1 for j e {I, . . ., p}, with Cj the jth standard unit 
vector in W\ can be imposed as follows. Given continuous functions Ai, . . .,Ap : 
Ap_i —> M verifying Aj(ek) = 6jk (Kronecker delta) for all j, fc e {l,...,p}, define 



1 



1 ^ (I 



- 1 



Al^iw) A^(w) tf ^ {A]:{ej) 
logAjf (iv) = log Af ^Civ) - ^ Aj{w) logAf °(e,). 



(3.1) 
(3.2) 



;=i 



In case of known margins, variance-minimi zing weight functions Aj can be deter- 
mined adaptively by ordinary least squares (|Segersl . 120071 : iGudendorf and Segers . 
201 ih . However, if the marginal distributions are unknown, these endpoint cor- 



rections are asymptotically irrelevant (IGenest and Segersl . 120091) . since 

iogA-(.,)=-i|:iog(iog(^)).];'iog(iog(i))d. 



o 



{\ognf 



as n ^ oo. Nevertheless, in finite samples, the simple choice Aj{w) = Wj can make 
quite a diff"erenc e. Similarly, for unknow n margins, the multivariate extension of 
the estimator in iHall and Tajvidil (|2000h simplifies to A^'^(w) = Al{w)/Al{ej) = 
A»{l+0(n-'logn)}. 

The next lemma establishes a functional relationship between A^ and A^^*^ on 
the one hand and the empirical copula C„ on the other hand. Recall the empirical 
copula process C„ in (12.71) . 



Lemma 1. For w e A„_i, we have 



u 



jii/2 / 1 3_\ ^ f £ (j^m 

\Al(w) A(w)l Jo "' 

ni/2(iogA,fV)-logA(H'))= r C„(""',...,- . , 

Jo u log u 



,m"'")- 



dw 



(3.3) 
(3.4) 



The proof is not diffe rent from the one in dimension two and can be found in 
Genest and Segersl (|2009l . Lemma 3.1). Equations (13.31) and (13.41 ) are instrumental 
for proving the weak convergence of the processes 



Al = n^'HAl-A) and 
Theorem 1. Under AssumptionUlabove, Aj, 



CFG 



,l/2/?CFG 



^i/Z(-^^ru _ ^y 



and . 



,CFG 



,CFG 



with 



AV) = -A^(w) C(m^\...,m"''') — 
Jo w 

C(m^'',...,m">) / , 

M lOgM 

as n ^ oo in the space ^(Ap_i) equipped with the topology of uniform conver- 
gence. 



9 



The main idea of the proof consists in substituting C„ in (13.31 ) and (13.41 ) by 
Stute's expansion and to conclude using a refined version of t he continuous map 



ping th eorem. As the proof follows the same lines as the one in lGenest and Segers 



(|2009|) . we will just point out the main adjustments. 



4. Projection estimator 

The estimators of the Pickands dependence functions considered so far are in 
general not valid Pickands dependence functions the mselves. In this section, we 



adapt the methodology in iFils-Villetard et al.l (|2008h to project a pilot estimate 



An onto the set s^ of Pickands dependence functions of p-variate extreme-value 
copulas. To this end, we view ^ as a closed and convex subset of the real Hilbert 
space L^(Ap_i) with Ap_i equipped with (p - l)-dimensional Lebesgue measure 
when viewed as a subset of M'' ' . The inner product and the norm on L^(Ap_i) are 
denoted by {f,g) = f fg and Wfh = i{f,f)y'^ respectively. 

The orthogonal projection of an initial estimator A„ for A, for example the 
Pickands or the CFG estimator, onto £/ is then defined as 

AP' = U{An\£/) = argminllA - AJI2. 

Projections being contractions, it follows that HAP'^ - AII2 ^ ||A„ - AII2 for all A G 
£/: the L^-risk of the projection estimator is bounded by the one of the initial 
estimator. 

For practical computations, we are obliged to refer to finite-dimensional sub- 
classes £/,„ c £/, yielding the approximate projection estimator 



fpi- 



Ii{AnWm) = argminllA -A„||2. (4.1) 



Aes/„, 



For each positive integer m, the class s/m will be defined as the set of Pickands 
dependence functions characterized by discrete spectral measures H with fixed 
and finite support depending on m. 

Specifically, let "Vp^m be the (finite) set of points v = (vi, . . . , Vp) e A;,_i such 
that kj = mvj is integer for every j e { 1 , . . . , p}, so that in fact v = (k\ /m, . . . , kp/m) 

where kj e {0, . . . , m} and ki -\ \- kp = m. The cardinality of Yp^m is of the order 

0{mP~^) as m ^ 00. Let M'p be the set of spectral measures on Ap_i and let 
^p^m be the set of (discrete) spectral measures H e J^p supported on 1^,,,,, that 
is, H = Jjy^r K <5v, with 6^, the Dirac measure at v and where h^ = H({v}) is the 



' p,m 



10 



spectral mass of the atom v. The vector h = (/Zv)ver„„ satisfies the constraints 

{ (4.2) 

[ Y.ver,,„,Kvj= 1, Vj €{!,...,/?}, 

the second constraint stemming from (11.31) . 

The Pickands dependence function A of a spectral measure H in J^,„j can be 
written as 

Ah{w) = ^ hy max{wivi, . . .,WpVp}, w e Ap_i. (4.3) 

Being a linear combination of piecewise linear functions, the function A in (14.31) 
is itself piecewise linear. All Pickands dependence function of the form (14.31) 
will be collecte d in the class Mn- The ne xt result can be seen as the equivalent 



of Lemma 2 in lFils-Villetard et al.l (|2008h stating the denseness of the piecewise 



linear Pickands dependence functions. 

Lemma 2. For every H e J^p and every positive integer m, there exists H^ £ 
J^p^in such that the Pickands dependence functions A and A,„ ofH and H„, respec- 
tively satisfy 

sup \A,n{w)-A{w)\^^. (4.4) 



The bound in ( 14.41 ) implies that sup^^^ inf^e<, 11^ ~ M\i = 0(m ^) as m ^ oo 



This r ate is perhaps not sharp, for in case p = 2, Lemma 2 in lFils-Villetard et al 



(l2008h states the rate 0(m~^'^). It remains an open problem whether the latter rate 
can also be achieved in general dimension p. 

la practice, the task is to compute the vector h such that the function 

A'^^iw) = Aj^iw) = ^ hy maxfwivi, . . . , WpVp}, w e Ap_i, 



solves (14.11) . The vector h is given as the solution to the least-squares problem 

h = argminllAft - A„||2 = argmin(||Aft||2 - 2(Aa,A„)), (4.5) 

h h 



11 



with h subject to the constraints (I4.2I ). The optimisation problem in (14.51) is a 
quadratic program with linear constraints, which in matrix notation reads 



cTh I, subject to 



/I ^ 
h = arg mini -h Dh 
h ^2 



The matrix D and the vector d regroup all the scalar products of the form 




(4.6) 



/ 

J^„ 



max(w v) max(H' v') &w and 



/ 



max(H'v)A„(H')dH' 



respectively, for v,v' e Yp^m. The p equality constraints 2ver,,„^vV; 
encoded by means of the matrix C and the vector c. 



1 are 



Fo r implementation, we used the R-package quadprog (|Turlach and Weingesseli 
20101) for solving quadratic programs under linear constraints. Although there ex- 
ist multiple packages for numerical multivariate integration, we preferred to com- 
pute all the integrals appearing in D and d using Riemann sums on the same fine 
grid. By so doing we reduce the risk of numerical problems as we impose D to be 
positive definite. 

The de rivation of the asyrnptotics of the projection estimator follows the same 
lines as in Ipils-Villetard et all pOOSb . Assume that e,;\A„ - A) -^ ^ in L^(Ap-i) 
where ^ is a random process in L^(Ap_i) and < e„ ^ 0; for the Pickands 
and CFG estimators, we have e„ = n~^^^ and we have weak convergence with 
respect to the uniform topology, which implies convergence with respect to the 
L^-topology. 



By Lemma 1 in lFils-Villetard et al.l (|2008|) . we have 



\\A!:-An\2^[6,„{2\\A„-An\2 + 6J] 



1/2 



with Sm = IIAP'' - n(AP''|^M)ll2- From Lemma |2] above, we have 6,n = 0(1 /m) as 
m ^ oo. As a consequence, i f m = m„ is such that 1 /m „ = o(s„) as n ^ oo, then 
ll^m -^P^b = op(e„). From iFils-Villetard et all (I2008L Theorem 1), we deduce 
that 



s;,\A^: - A) = s-„\A^' - A) + op(l) ^ n(^|T^(A)) (n ^ oo) 



(4.7) 



in the space L^(Ap_i), where T_^(A) is the tangent cone of £/ at A, defined as the 
L^-closure of {A(A - A) : A ^ 0,A e £/}. 

Interestingly, equation (14.71) implies that the choice of m is not to be seen as 
a bias-variance trade-off problem but rather as a discretization problem. As soon 



12 



asm = m„ converges to infinity faster than !/£„, the finite-dimensional projection 
estimator AjJ[ has the same limit behaviour as the 'ideal' projection estimator AP''. 
In practice, we will choose m sufficiently large so that any further increase of 
m does not make any significant difference, of course subject to constraints on 
computing time and numerical stability. 

Finally, note that the convergence in (14.71) is with respect to the L^ -topology 
only, even if originally the weak convergence of s~^(A„ - A) took place in the 
stronger ^""-topology. The asymptotic distribution of the projection estimator un- 
der the ^""-topology remains an open problem. 



5. Simulation study 

A simulation experiment was conducted to compare the finite-sample perfor- 
mance of the following four estimators: 
PD - the end point-corrected Pi ckands estimator in (13.11) with Aj(w) = Wj, in the 

spirit of lOeheuvelsl (|l99lh : 
PD-pr - the projection estimator in (14.11 ) with the previous end-point corrected 

Pickands estimator as initial estimator; 
CFG - the endpoint-corrected CFG estimator in (13.21 ) with Aj(w) = wf, 
CFG-pr - the projection estimator in (14.11) with the previous end-point corrected 

CFG estimator as initial estimator. 



The set-up of the expe r iment was as follows. Following IZhang et al.l (|2008|) 
and lGudendorf and SegersI (|201 lb . random samples were generated from a trivari- 
ate extreme- value distribution with asymmetric logistic dependence function A 
(lTawnLll99(3h : 



A(w) = (e^'^w]'" + (f> 



l/o- 1/q'\Q' 



ll/ff,^ l/ff 



l/a 1/Q'\a 



^w;'T + (6''' < +0""W3 



r 



I /a 



+ {o'l^w''" + (p'i"w\'"f + ijj{w\'" + w;'" + w''"f + \-e-(p-ijj, (5.1) 



for w e A2, with parameter vector (a, 9, (p, ijj) e (0, 1] x [0, 1]^. For this model. 
Assumption [2] can be verified by direct calculation. The dependence parameter 
a ranged from 0.3 (high dependence) to 1 (independence, A = 1) and the vector 
{(p, ip, 6) was set equal to either (0, 1 , 0) (symmetric logistic copula or Gumbel cop- 
ula) and (0.3,0,0.6) (an asymmetric logistic copula). For each distribution, 1000 
samples were generated of " size n e {50, 100 , 200}. Simulations were performed 
using the R- package evd (IStephensonl |2002|) . which implements the algorithms 
presented in IStephensonI (|2003h . The discretization parameter m was set to 20, at 
which value the grid ^^20 contains 231 points. 



13 



Monte-Carlo approximations for the mean integrated squared error (MISE) 
E[ J(A - A)^] for the four estimators considered above are reported in the tables 
below. The three main findings are the following: 

1- The projection step yields a gain in efficiency, especially in case of weak de- 
pendence. 

2- Without projection step, the CFG estimator outperforms the PD estimator. 

3- After the projection step, the PD-pr estimator is more efficient than the CFG-pr 
estimator in case of independence and weak dependence (a ^ 0.9), but less 
efficient otherwise (a ^ 0.7). 

Further, as the dependence increases, all estimators tend to perform better. In 
accordance with asymptotic theory, the MISE is roughly proportional to 1 /n. 
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Appendix A. Proofs 

Appendix A.l. Proof of Proposition\l\ 
If M e (0, 1]^, then - log m e Wp,y and 

C{u) . 
Cj{u) = ^-^tj{-\ogu). (A.l) 

Uj 

The assumptions on € imply continuity of Cj on the set (0, 1]''. If m e [0, \y with Uj > 
and Ui = for some / € {I, . . .,p] \ (j), then Cj{u) = and continuity of Cj at such u 
follows from the fact that ^ ^y ^ 1 and ^ C(v) ^ min(v). 
If (L2) holds, then also 

sup max(x;, Xj) \{ijix)\ < cxj, 

xeW,,jnWpj 

that is, without the condition x\ + ■ ■ ■ + Xp - 1. This result is based on the fact that the 
function £ is homogeneous of order one: ^(^jc) = st{x) for all s € [0, oo) and x e [0, oo)P . 
Hence if tj exists on Wp_j, then for all s € (0, co) and x £ Wpj we have 

d d 

s€j{sx) = - — £{sx) = - — sf{x) = sij{x) 

OXj OXj 

and thus 

ij(sx) = ijix). 

Taking partial derivatives again, we find for all s € (0, oo) and x e Wpj n Wpj that 

sfijisx) = -—ij{sx) = -T-ijix) ^ 'iijix) 

OXj OXj 



and thus 



'ijiisx) - s-^ 'iijix). 



It follows that 

max(5;c;, sxj) 'iij(sx) = max(jc/, xj) ijjix), 
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that is, the map x i-> max(x;, xj) 'tij{x) is constant on rays through the origin. 

Next, we show the equivalence of (L2) and (C2). Fix /, 7 e { 1, . . . , p}, not necessarily 
distinct and let u e Vpj n Vpj. On the one hand, if m e Vpj n Vpj n (0, 1]'^ (meaning that 
every component of u is different from 0), then - log u e Wpj n Wpj and 



Cu(u) 



UiUj 
Uj 



with the convention that the partial derivatives of i are evaluated in - log m. On the other 
hand, if a e iVpjr\Vpj)\{0, ly (i.e. at least one coordinate of m vanishes), then Cij(u) = 0. 

We have to verify two things: first, the continuity of C/y at points in the set {Vpj n 
Vp.j) \ (0, 1]''; secondly, the finiteness of the supremum in (C2). 

First, let u e Vpj n Vpj n (0, 1]''. Let ^ be a positive constant not smaller than the 
supremum in (L2). By assumption (L2) and the fact that ^ /'y ^ 1, we have 

min(H) / K 



|C,(„)K]I^H1^(i, 

U;U i \ 



max(- log Uj, - log My) 



Continuity of C,y at points in the set Vpj n Vpj \ (0, 1]'' follows. 

Secondly, as min(M)/(M,My) ^ min(l/M,-, l/uj) and - log x ^ 1 - x for all positive x, 

|C,/(«)I ^ min(— , —]{l+K min( 



\ Uj UjJ \ \ 1 - M; ' 1 

/I 1 \ 



/I 1 \ /I 1 

{\+K) mm — , — mm , 

\ Uj Uj/ \l — Uj 1 - Uj 

( ' ' \ 

\Ui{\ - Ui) Uj{\ - Uj)/ 



^ (1 + ^) min( , ) 

\U,(\ - Ui) U,i\ - U;)/ 



which is equivalent to condition (C2). 

Appendix A.2. Proof of TheoremU} 

The proof of theorem [J will require the following preliminary result on weighted 

empirical copula processes. Recall the process a,j in ( 12.21 ). Define qg(t) = t^(l - t)^ for 

t e (0, 1) and a fixed value 6 € (0, 1/2). Write E = (0, 1]'^ \ {(1, . . . , 1)). Define the process 

Gn,0 on [0, \y by 

{ ""^") ,, ifweE, 
G„e(«)-r'^™"^"^^ (A.2) 

\0 ifMe[0,lF\E. 

Similarly, define the process Ge on [0, 1]'' by replacing a„ in (IA.2I) by its weak limit a. 



see (12.101) . The following result generalizes Theorem G.l in lOenest and SegersI (|2009|). 
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Lemma 3. For every 6 e (0, 1/2), the trajectories of Go are continuous almost surely and 

Gn,e -^ Ge in r{[0, lY). 



Proof (Lemma |3]|. The proof is entirely analogue as the one of Theorem G. 1 in lGenest and Segers 



(120091) ■ For completeness, we sketch the main lines. 
Fix M € E and define the mapping /„ : E — > M by 

, , , l(0,«](s) - C{u) 

fu{s) - — — ■ , s e E, 

^6i(mm(M)) 

and consider the class 

^ = {/„ : H e E) U {0}, 

where of course stands for the zero function. The space ^ will be endowed with the 
metric 

p\f,g) = 'P{f-gf f,g^^. (A.3) 



Here, we adopt the notations of Ivan der Vaart and Wellnerl (119961) : P denotes the proba 



bility distribution on E corresponding to C and Pn denotes the empirical measure of the 
sample (Un,. . . , Ujp) for / e { 1 , . . . , ?i), that is 

/I " 
fdC, ¥nf = -yKUu,-..,Ui,p). 

Moreover, put G„ = ?i'''^(Pn - P), viewed as a random function on ^. 

We will show that the collection ^ is a P-Donsker class, i.e. there exists a P-Brownian 
bridge G such that 

G„ -^ G in ^~(^) as n ^ oo. 



It is s ufficient to verify the conditions of Theorem 2.6.14 of Ivan der Vaart and WeUner 



(Il996l) . The function f on E defined by 



F(su. ..,sp) = pmax|^7^ . . . , s;f,{l - siT^ ... ,(1 - SpT^^} 
is a suitable envelope function for ^. The fact th at ^ is a VC-major class an d is pointwise 



sepai^able follows from the same arguments as in lGenest and SegersI (120091) . 

For the moment G„ is defined on ^ with the metric p in (IA.3! ). Consider the map 
4> : [0, 1]^ -^ ^ defined by ^(m) = /„ if m £ E and ^(m) = if m e [0, 1]^ \ E. Then 
Gnfi = G„o0 and Ge = Go0. The map £°°{^) -^ ^°°([0, 1]'') : z h^ zo,^ being continuous, 
the continuous mapping theorem permits to conclude that Gnfi '^ Go in ^°°([0, 1]''). Since 
the trajectories of G are p-continuous almost surely and since is continuous, it follows 
that the sample paths of Go are continuous almost surely as well. This concludes the proof 
of Lemma [3] 
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We now proceed with the proof of Theorem [T] Define 



for w e Ap_i. Applying the change of variables u = e^* in Lemma [1] we find that the 
processes Bf^ and B^^'^ can be written as 

Bn(w)= C„(e-"'^...,^-"'''0/^Wd5, (A.4) 



in terms of a function h on (0, oo) which is h^is) = 1 for the Pickands estimator and 
h^^^{s) = l/s for the CFG estimator. In what follows, the function h denotes either h^ or 

Put /„ = l/(n + 1) and k„ = plog{n + 1) and split the integral on the right-hand side 
of (IA.4I ) into three parts: 

Hn rK n°° 
Bn{w)= + + = /l,„(H') + /2,„(H') + /3,„(H'). (A.5) 



n *^ '^n 



We will first prove that with probability one, the first and the third term on the right-hand 
side converge to zero uniformly in w. 

• \f s e [0, In], then e'' ^ 1 - /„ and thus C„(e""'", . . . , e^'^'p") = 1, which implies 

^ h^„{w) = f " n^'^il - e-^^(»')) h{s) ds 
Jo 

^ n^'^ [ sh{s)ds ^ n^'^ln ^ n-^'^. 
Jo 

• If 5 e {kn,oo), then Wj ^ l/p and thus e~^^'J'^ < l/{n + 1) for at least one j e 
{1, . . . ,p}, so that C„(e""'", . . . ,e~"'"') ^ 0, which imphes 



f 

Jk., 



lk„ 

„l/2 



A{w) 
where we used the fact that A{w) ^ max(H') ^ l/p. 
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As a consequenc e, the only non -negligible term in (IA.5I ) is h^n- By Assumption [T] and 
by Proposition 4.2 in lSegersI (120111) . Stute's expansion (I2.8l) -( |2!9l ) is valid, so that we can 



write 

p 

j=i 

where 

/0,n(H') - r " «„(£-"' ^ . . . , £>-*"■') h(s) ds, 

JjAw) - r " anjie-'^n Cj{e-^'\ ..., ^-">0 Ks) As, je{\,...,p}, 

Jl„ 
rk„ 
Jp+U^) = R„ie-''' ',..., e-">0 his) ds. 

In view of the bound ( I2.9I ) on R„, the term 7p+i,„ is negligible: as n —> oo, 

sup |/„+i „(w)\ = 0{n~^^^ login) f " his)ds) -^ 0, almost surely. 

Fix 6 e (0, 1/2) and recall the process G,,,^ in ( IA.2I ). We have 

/an(H') = Gn,eie~''"' ,■■■, e-'''-")Kois, w) his) ds, 

Jin 

Jj,niw) = f " G„,e(l, . . . , l,e-•™^ 1, . . . , l)Kjis,w)his)ds, j e {1, . . .,p], 

Jin 

with 

Kois, w) - 9e(min(e-™'' , . . . , e"™")), 

i^:/^, iv) - qeie-'''^ ) C/e-™' , . . . , e-^"'0, je{l,...,p}. 

The functions Z^o^ ■ ■ ■ ,Kp satisfy the bounds 

^ Kjis,w) ^ Kis) - /l(o,i](5) + e-^^'P^' l(i,co)(^), 

;• e {0, . . . , p), 5 e (0, oo), w e Ap_i . (A.6) 

To prove these bounds, use equation (|A. IK the fact that ^ ^y ^ 1 and ^ C(v) ^ 
min(v) for v e [0, 1]^ and the fact that max(M') ^ 1/p for w € Ap_i. The function K in 
(IA.6I ) satisfies j^ Kis) his) ds < oo. 
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By Lemma[3]and the extended continuous mapping theorem (Ivan der Vaart and Wellned. 

19961 . Theorem 1.11.1), we find 



7=1 

in ^"(Ap_i), where 

r*00 

Joiw) = Ge(e--™'' , . . . , e"™'') ^o(^, h') h{s) ds, 
Jo 

Jj{w)= Geih...,l,e-'"'\\,...,l)Kjis,w)h(s)ds, je{l,...,p}. 
Jo 

Substituting the definitions of the process Go and the functions Kq, . . . , Kp, we obtain 



-f 

Jo 



B{w)= C(e~"'' ',..., e'">0/jWd5, weAp^i. 
Jo 

An application of the functional delta method (Ivan der Vaart and Wellnen. 1 19961 Theo- 
rem 3.9.4) now yields the result. 

Appendix A.3. Proof of Lemma^ 

The proof is consti'uctive and consists of the following steps: 

1. Construction of the spectral measure H,„. 

(a) Discretisation of H yielding a measure G,„ on Ap^i, which is not necessarily a 
spectral measure. 

(b) Modification of G,„ into a genuine spectral measure H,„ . 

2. Proof of the inequality (l44l ). 



1. Construction of tlie spectral measure Hm- For v € 1^,,„> consider the set Ap_i v,m of 
points t e Ap_i such that Vj ^ tj < vj + l/m for every j e {I, . . . , p - l); recall that 
tp = I - ti - ■■■ - tp-i, so that necessarily Vp - (p - \)/m < fp ^ Vp. The collection of 
sets {Ap_i v,m : v € ^,,m) constitutes a partition of Ay,_i. Indeed, for every point t e Ap_i 
there is a unique point v e 1^,„, such that t e Ap_i^v,m: Let vy be the integer part of mtj for 
je{l,...,p- l}andputvp ^ 1 -vi Vp_i. 

(a) Discretisation of H, yielding G,„. Define a discrete measure G,n on Ap_i with 
support contained in ^,,„ by Gmiiv]) = //(Ap_i ,,,„) for v € Yp^m- In words, the mass 
assigned by the spectral measure H on the set Ap_i v „, is relocated to the corner point v. 

Since the sets Ap_i v„, constitute a partition of Ap_i, the total mass of Gm is still 
Gm(Ap_i) = H(Ap-i) = p. However, Gm does not need to verify the moment constraints. 
For 7 € { 1, ...,;?} we have 

j tjdG,n{t)= Yj ^j(^'ni{v})= 2 Vy//(Ap^i,v,,„), 



Ap-l 
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which in general is not equal to unity. 

Still, the moment constraints are not far from being verified. For t e Ap_i v„, and 
j e {\, . . . ,p - 1 ) we have Vj ^ tj < Vy + 1 jm. Integrating these inequalities over 
t e Ap_i „ „, with respect to H and summing them over v G 1^ ,„ yields 

f tj dCnit) < 1 < r tj dG,n{t) + -, je{l,...,p-l}. 

As a consequence, there exist numbers Cj e [0, 1) such that 



I 



tjdG,„it) = l-^-^, je{l,...,p-\]. 



A , m 



(b) Modification ofGm into a spectral measure H,n. We will modify G,„ into a genuine 
spectral measure H,„ by (slightly) increasing the masses at the vertices e\,..., Cp-i, where 
ej is the 7th coordinate vector in M''. Specifically, we set 



//„, = (1 - ao) G„, + ai6ei+--- + ap-x 5, 



p-i 



for some nonnegative numbers a^, . . . ,ap^i to be determined by the moment constraints. 
For 7 e { 1 , . . . , p - 1 }, we must have 

1 = I tjdH,„{t) = (1 - aQ){\ - Cj/m) + aj. 
Ja,,_, 

In addition, the total mass must be equal to 

p ^ H,„{Ap-i) = (1 -aQ)p + ai + ■■■ +ap^i. 

Substituting a^ = 1 - (1 - ao)(l - Cj/m) into this equation and solving for ao yields, after 
some algebra. 



flo 



W + ZfJi Ci 
Cj + Zfj/ Ci 



m + Y.U 



Ci 



ie{l,...,p-\]. 



This concludes the construction of the spectral measure Hm- Note that ^ a^ < pim for 
every; £ {0,...,p-\}. 



2. Proof of the inequality (14.41 ). For w,t e Ap^i, write 

f{w,t) ^ max{wi?i,..., Wpfp}. 
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The Pickands dependence function A,„ of the spectral measure H,„ constructed above is 
given by 



Aniiw) 



f{W,t)H,nidt) 



-L 

^(I-uq) I f{w,t)G,n{dt) + aiWi + --- + ap-iWp-i. 

JA,,-i 

Put Bm{w) = //(h', t)Gini<it), the "Pickands transform" of G,„. Clearly Bm ^ and 
B,„ is convex, being a weighted average (over t) of the convex functions w i-> f(w, t). 
As a consequence, B,„{w) ^ max{B„,(^i), - - - , B,„{ep)]. Now B,„{ej) - jtj Gmidt), which 
is equal to 1 - cj/m if j e {I, . . . ,p - 1} and which is equal to p - Tjf=i (1 ~ c,/m) = 
1 + l^fr/ a/m = 1/(1 - ao) if j = p. It follows that B„^{w) ^ 1/(1 - ao) for all w e Ap_i. 

We obtain, on the one hand, 

A,„{w) ^ B,„{w) + aiwi + ■■■ + ap-iwp-i 



^ Bm{w) + max(ai, . . . , Up-i) < B,„{w) + — 



m 



and, on the other hand. 



1 P-i 

Uq 1 v~i 

Am{w) ^ (1 - ao)B,n{w) ^ B,„{w) - = B,n{w) > : 



1 -ao 



m 

>Bm{w)-^. 

m 



i=\ 



Therefore, 



\A{w)-Al{w)\ ^ \A{w) - B,n{w)\ + \B,„{w)-A,n{w)\ < \A{w) - B^{w)\ + ^. 



m 



Furthermore, 

\A{w) - B^(w)\ ^ 



f{w,t)G,n(dt) 



Y^ f f{w,t)H{dt)- f 

'= ^ pju ' ' 

= 2 r f{wJ)H{At)-f{w,v)H{Kp^,^m) 

" Z r {f(}v,t)-f{w,v))H{At) 

2 r \f{w,t)-f{w,v)\H{dt). 
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By checking all possible cases one verifies that | max{a\ , 0:2) - ma.x(b\ , b2)\ ^ maxdai - 
Z71I, 1(32 - ^2!) for all real ^1,^2, b\,b2- An induction argument yields | max(c[i, . . . ,ak) - 
rmx{bi, . . . ,bk)\ ^ max(|ai - b\\, . . . ,\ai^ - b^\). It follows that 1/(^,0 - /(w,v)| ^ 
maxd^i - vil, . . . ,\tp - Vp\). As a consequence, 

\A{w) - Bm{w)\ < > I ^^ H{dt) = ^'' ~ ^^^ 



Inequality ( I4.4I ) follows 






24 







a = 0.3 


a = 0.5 


a = 0.7 


Q' = 0.9 


a= 1 


n = 50 


PD 


1.40-10-4 


5.44 • 10-4 


1.36 • 10-^ 


2.68 • 10-3 


3.44 • 10-3 




PD-pr 


1.37 


10-4 


5.14 


10-4 


1.21 


10-3 


2.08 


10-3 


2.44 


10-3 




CFG 


9.77 


10-5 


4.27 


10-4 


1.26 


10-3 


2.54 


10-3 


3.48 


10-3 




CFG-pr 


9.69 


10-5 


4.22 


10-4 


1.22 


10-3 


2.43 


10-3 


3.31 


10-3 


n = 100 


PD 


7.08 


10-5 


2.84 


10-4 


7.06 


10-4 


1.34 


10-3 


1.69 


10-3 




PD-pr 


6.99 


10-5 


2.74 


10-4 


6.53 


10-4 


1.03 


10-3 


1.08 


10-3 




CFG 


5.03 


10-5 


2.39 


10-4 


6.56 


10-4 


1.23 


10-3 


1.48 


10-3 




CFG-pr 


5.01 


10-5 


2.37 


10-4 


6.47 


10-4 


1.18 


10-3 


1.36 


10-3 


«-200 


PD 


3.31 


10-5 


1.43 


10-4 


3.92 


10-4 


7.02 


10-4 


8.71 


10-4 




PD-pr 


3.29 


10-5 


1.40 


10-4 


3.73 


10-4 


5.72 


10-4 


5.14 


10-4 




CFG 


2.45 


10-5 


1.23 


10-4 


3.39 


10-4 


6.40 


10-4 


6.56 


10-4 




CFG-pr 


2.45 


10-5 


1.23 


10-4 


3.36 


10-4 


6.19 


10-4 


5.78 


10-4 



Table 1: Symmetric logistic dependence function, (0, ijj, 6) - (0, 1,0): Monte-Carlo approximation 
of the MISE of four estimators of A based on 1000 random samples 
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»-0.3 


tt-0.5 


Q'-0.7 


a = 0.9 


n = 50 


PD 


1.42 • 10"^ 


1.72 • 10-3 


2.20 • 10-3 


2.88 ■ 10-3 




PD-pr 


1.22 


10-3 


1.45 


10-3 


1.74-10-3 


2.01 


10-3 




CFG 


1.15 


10-3 


1.41 


10-3 


1.84- 10-3 


2.77 


10-3 




CFG-pr 


1.10 


10-3 


1.35 


10-3 


1.75 


10-3 


2.60 


10-3 


n = \00 


PD 


7.67 


10-4 


8.76 


10-4 


1.10 


10-3 


1.51 


10-3 




PD-pr 


6.77 


10-4 


7.70 


10-4 


9.00 


10-4 


1.08 


10-3 




CFG 


5.90 


10-4 


7.06 


10-4 


9.05 


10-4 


1.20 


10-3 




CFG-pr 


5.70 


10-4 


6.85 


10-4 


8.68 


10-4 


1.10 


10-3 


n = 200 


PD 


3.92 


10-4 


4.72 


10-4 


5.84 


10-4 


7.60 


10-4 




PD-pr 


3.52 


10-4 


4.30 


10-4 


5.08 


10-4 


5.19 


10-4 




CFG 


3.01 


10-4 


3.31 


10-4 


4.43 


10-4 


5.81 


10-4 




CFG-pr 


2.92 


10-4 


3.22 


10-4 


4.29 


10-4 


5.36 


10-4 



Table 2: Asymmetric logistic dependence function, {(p, i/r, ff) - (0.3, 1, 0.6): Monte-Carlo approxi- 
mation of the MISE of four estimators of A based on 1000 random samples 



26 



